Appearance of neutronization peak and decaying supernova neutrinos 



O 

o 

(N 



(N 

(N 
> 
O 

o 

in 
o 

O 

a: 
i 



X 



Shin'ichiro AndcQ 
Department of Physics, School of Science, The University of Tokyo, 
7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 
(Dated: Submitted April 13, 2004; accepted May 20, 2004) 

Nonradiative neutrino decay, which is not satisfactorily constrained, possibly and significantly 
changes the detected neutrino signal from galactic supernova explosions. We focus on the appearance 
of a sharp peak due to a neutronization burst in the time profile; this phenomenon would occur 
if the original v £ , produced at the neutrinosphere and becoming V2 or v-j, at the stellar surface, 
decays into a lighter antineutrino state such as v\ or vi coupled to v a . This is due to the fact that 
the signature of the neutronization burst is common to all numerical simulations, contrary to the 
spectral energy distribution of each flavor neutrino and antineutrino, which is still under intense 
debate. Therefore, the appearance of a neutronization peak in the P e signal, if it were detected, would 
clearly indicate the nonstandard properties of neutrinos; the nonradiative neutrino decay would be 
one of the possible candidates. Using a newly developed formulation that includes flavor conversions 
inside the supernova envelope and neutrino decay during propagation in a vacuum, we calculate the 
expected neutrino signal at the detectors; the lifetimes of three modes m, rrs, and T2S are taken to 
be free parameters. We further introduce simple quantities, which represent a peak sharpness of the 
time profile and spectral hardness, and investigate the parameter dependence of these quantities. 
As the result, it is found that they are quite dependent on the relevant parameters, but it would be 
quite difficult to distinguish models using the signal obtained by the Super-Kamiokande detector; 
the future megaton-class detectors would have sufficient sensitivity. We also compare the neutrino 
decay model with another mechanism — i.e., resonant spin-flavor conversion — which also may give 
the appearance of a neutronization peak, and conclude that these two independent mechanisms give 
a very different signal and one can be distinguishable from the other. 

PACS numbers: 95.85.Ry, 13.35.Hb, 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

In recent years, we have made great progress con- 
cerning our knowledge of neutrino properties; especially, 
many ground-based experiments, which observed atmo- 
spheric [lj, solar 0, and reactor neutrinos ||, have re- 
vealed nonzero neutrino masses and mixing angles — i.e., 
properties beyond the standard model of particle physics. 
However, many other neutrino properties are left un- 
known, such as the nonzero magnetic moment and neu- 
trino decay. Fortunately, our current knowledge of the 
mass differences as well as mixing angles enables us to 
consider these further exotic properties. 

The most stringent and precise limits on both the neu- 
trino magnetic moment and lifetime of nonradiative neu- 
trino decay are obtained by solar neutrino observations. 
The basic technique for both cases is as follows. We al- 
ready know that the famous solar neutrino problem is 
best explained by the large mixing angle (LMA) solu- 
tion. If other mechanisms, such as the magnetic moment 
or nonradiative decay, work in nature, it should change 
the observed signal slightly; i.e., no positive claim can 
then be used to put a limit on the other exotic prop- 
erties. For example, the neutrino magnetic moment, if 
ever, would produce T> e via spin-flavor precession inside 
the Sun 0,13 • However, the recent KamLAND experi- 
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ment report that no v e candidates were found from the 
Sun and this result is used to obtain the upper limit 
IU,„Bt(0-05Rq) < 10 _5 ^s G, where \ib is the Bohr mag- 
neton and Bt represents the transverse component of 
the solar magnetic field @ . As for nonradiative neutrino 
decay, the lower limit to the lifetime is obtained to be 
T12/1T1 > 10~ 4 s/eV, owing to no positive signature of 
the decay in the flux and spectrum of the solar neutrinos 
(see also Ref . 6] , for a more recent and stringent limit 
using the KamLAND data). For other laboraory bounds, 
we refer the reader to Refs. jaQ. 

These current limits are, however, still rather weak and 
such exotic mechanisms we consider potentially work in 
a more extreme environment as actually expected in the 
case of core collapse supernovae. Spin- flavor conversion 
of supernova neutrinos, induced by the interaction be- 
tween the neutrino magnetic moment and the supernova 
magnetic fields, has been studied by many researchers 
0, [n], 0, 0| and found to give a leading effect on the 
neutrino spectrum and luminosity curve detected at the 
Earth. Also, in the case of neutrino decay, because galac- 
tic supernovae are expected to be located at a much more 
distant place compared with the Sun, a far more stringent 
lower limit to the lifetime is expected. In addition, it has 
recently been proposed that the diffuse background of 
neutrinos emitted by past supernova explosions may be 
used to probe the lifetime of neutrinos |13|. Thus, core- 
collapse supernovae are considered to be a wonderful as- 
trophysical event, which can also be used as a laboratory 
for particle physics beyond the standard model. Other 
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high-energy astrophysical objects are also expected to be 
available for this purpose |14j . 

However, the expected galactic supernova neutrino sig- 
nal including nonradiative decay has not been discussed 
precisely; only very rough discussions have been done. 
Therefore, in the present paper, we give comprehensively 
the expected supernova neutrino signal at the large vol- 
ume water Cerenkov detectors on the Earth using realis- 
tic models of original neutrino spectrum and luminosity 
curve numerically calculated by Thompson et al. 
In particular, we focus on the case in which the sharp 
peak of a neutronization burst appears in the detected 
v e signal at the water Cerenkov detectors. This is due to 
the original u e , produced by neutronization of the pro- 
toneutron star matter, possibly decays into v e component 
during its propagation. In fact, in the previous paper we 
have pointed out that spin-flavor conversions in the su- 
pernova envelope can also cause the same phenomenon — 
i.e., th e ap pearance of a neutronization peak in the v e 
signal |Ty| (see also Ref. |12|). Hence, we compare 
the expected neutronization peak due to neutrino de- 
cay with that due to spin-flavor conversion and discuss 
their difference. Anyway, the appearance of such a signa- 
ture clearly indicates a nonstandard neutrino property; 
in that case, neutrino decay would be one of the possible 
candidates. On the other hand, the obtained neutrino 
spectrum would be useful, but the shape of the original 
spectrum is still matter of debate. Although the differ- 
ence among several groups is not very large, such a slight 
difference gives a large uncertainty when we discuss fla- 
vor conversion mechanisms or decay during propagation; 
i.e., whether the observed signature comes from the in- 
trinsic or extrinsic effect (such as decay) would be quite 
unknown at present. 

It should be noted that in this study we consider only 
vacuum neutrino decays. It is possible, however, to con- 
struct models where fast invisible decays can be triggered 
by matter effects [S E E E3 at the very high den- 
sities characterizing the supernova neutrinosphere, even 
in the absence of vacuum decays. In such scenarios, 
matter-induced decays (or interactions) might thus oc- 
cur before flavor conversions in supernovae, leading to a 
phenomenology rather different from the one considered 
in this paper. We emphasize that the results discussed in 
the following sections are generally applicable to vacuum 
neutrino decay occurring after flavor transitions and our 
approach is not constrained from the supernova cooling 
discussion as detailed later. In addition, with the lifetime 
considered in this study, the flavor conversions occur well 
before the decay and mass eigenstates are expected to be- 
come incoherent. Thus, the interference effects between 
the two phenomena as discussed in Ref. 2Q] can be ne- 
glected in our treatment. 

This paper is organized as follows. In Sec. [HJ we give 
several descriptions of models of decaying neutrinos and 
introduce a specific formulation of the decay rate from 
the Lagrangian. In Sec. IIIII an adopted model for the 
original (which means before occurring extrinsic processes 



such as flavor conversions or decays during propagation) 
spectrum and luminosity curve of supernova neutrinos 
are introduced; the effect of flavor conversions (without 
decay) is also described briefly. A newly developed for- 
mulation including both flavor conversions and decay is 
presented in Sec. II VI and we further give the calculated 
results obtained with the formulation in Sec. [3 Finally, 
we discuss several other possibilities in Sec. lVII and a brief 
summary is given in Sec. IVIll 



II. MODELS OF DECAYING NEUTRINOS 

In this paper, we study nonradiative two-body de- 
cay of the "invisible mode" — i.e., decays into possibly 
detectable neutrinos plus truly invisible particles — e.g., 
light scalar or pseudoscalar bosons such as the Ma- 
joron 21] . On the other hand, radiative neutrino decay 

— > has been considered in many papers (see Ref. 
|22j and references therein) and very stringent limits on 
the lifetime-to-mass ratio, r/m > 10 20 s/eV, have al- 
ready been set by several arguments [2i|. The Majoron 
models that cause nonradiative neutrino decay typically 
have tree-level scalar or pseudoscalar couplings of the 
form 



C = gijViVjX + >>.ji'j'r,i'. \ + H.c, 



(1) 



where x represents a massless Majoron, which does not 
carry a definite lepton number. For the coupling specified 
by Eq. ©, the decay rates into neutrino and antineutrino 
daughters are given by 0, 0] 
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where x = 7712/ mi, and we have dropped the subscripts 
on the coupling constants. Analogous expressions hold 
for v decay, with the replacements v «-> v. The decay 
widths in this paper are defined in the laboratory frame, 
and therefore the relation to the rest-frame lifetimes is 
r(E) = m/(rE). From these two expressions (0) and 
© it is easily seen that the decay rate is dependent on 
whether the helicity flips or conserves as well as whether 
the neutrino masses are quasidegenerate (x m I) or hi- 
erarchical [x ^> 1). In the case of strongly hierarchical 
masses, we obtain r u2 ^ Ul « r^ 2 ^ Pl ; on the other hand, 
in the case of quasidegenerate masses, Eqs. © and © 
lead to the relation 7V 2 _» Vl 3> r^—t^, unless the cou- 
pling constant g is strongly suppressed compared with 
h. Therefore, if the neutrino masses are quasidegenerate, 
the produced neutrinos decay into neutrinos dominantly 
(helicity-conserved channel), while little into antineutri- 
nos (helicity- flipped channel); hierarchical masses allow 
both decay channels by the same branching ratio. 
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The distribution of the energies of daughter 
(anti)neutrinos is very dependent on whether the 
masses are hierarchical or quasidegenerate. In the for- 
mer case, the energy of daughter neutrinos is generally 
degraded, and its distribution is given by 

^ c (E p ,E d ) = ^, ^ f { Ep ,E d ) = ^(l-^J, (4) 

where E p and Ed represent the energy of parent and 
daughter neutrinos, respectively, and ij; c and ipf are 
the distribution function of the helicity-conserved and 
hclicity- flipped channels, respectively. In the latter case, 
on the other hand, the daughter neutrino energy is al- 
most the same as that of parent neutrinos; the energy 
distribution becomes a delta function S(E P — Ed). 

As already mentioned in Sec. HJ we are mainly inter- 
ested in the appearance of a neutronization peak at de- 
tectors, which dominantly catches P e ; this is because the 
signature of the neutronization burst is very common 
to any supernova simulations, and its detection in the 
9 e channel would be a smoking gun to the nonstandard 
properties of neutrinos. Thus, the case of quasidegener- 
ate masses, which prohibits the helicity-flipped channel — 
i.e., v e — > v e , is not our prime interest. Although the ob- 
tained spectrum would also be helpful even in that case, 
we assume that the neutrino masses are strongly hierar- 
chical (mi <C ?7i2 <C 777,3) from this point on. 

At the end of this section, we mention the constraints 
on the coupling constants g and h from the supernova 
cooling argument. In the dense core of supernovae, con- 
trary to the decay in vacuum, the Majoron cooling pro- 
ceeds mainly via v e v e — > x- A conservative upper limit 
on the coupling constant \g ee \ is obtained by the fact 
that the luminosity of the Majoron, L x , should not ex- 
ceed 3 x 10 53 erg s _1 0]> because we already know that 
almost all the gravitational binding energy of new-born 
neutron stars is released as neutrinos by the actual de- 
tection of supernova neutrinos by Kamiokande [2j| and 
1MB 2G] detectors. This discussion translates into the 
bound I p ee I <4x 10~ 7 and it is the strongest constraint 
on the parameter compared with other experimental ones 
such as of neutrinoless double-beta decay with Majoron 
emission [27j. Bounds on other parameters such as \g efJl \ 
and \gfifil are also obtained by the similar arguments (for 
a detailed discussion, see, e.g., Ref. 01)- Our discus- 
sion in the present paper, however, is completely free of 
such stringent constraints. This is because the relevant 
parameters in our case are the coupling constants in the 
basis of mass eigenstates gtj, while those in the cooling 
argument are in the basis of flavor eigenstates g a p- These 
two expressions in different bases are related to each other 
as 

9ij = 229o t pUl i Uf} :j . (5) 

a/3 

Since the mixing matrix U a i contains several unknown 
parameters such as #13 or the CP- violating phase 5, the 



limits from supernova cooling do not directly translate 
into gij relevant for our study. 

III. ORIGINAL SUPERNOVA NEUTRINO 
SIGNAL AND FLAVOR CONVERSIONS 

We adopt as the original neutrino spectrum as well as 
the luminosity curve, the results of the numerical simu- 
lation by Thompson et al. ^5|; we use the model calcu- 
lated for the 15Mq progenitor star. Their calculation has 
particularly focused on shock breakout and followed the 
dynamical evolution of the cores through collapse until 
the first 250 ms after core bounce. They have incorpo- 
rated all the relevant neutrino processes such as neutrino- 
nucleon scatterings with nucleon recoil as well as nucleon 
bremsstrahlung; these reactions have recently been rec- 
ognized to give a non-negligible contribution to the spec- 
tral formation. In Figs. ^ and [21 we show the original 
luminosity curve and number spectrum of neutrinos, re- 
spectively. In these figures, u x represents nonelectron 
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FIG. 1: Original luminosity of the emitted neutrinos as a 
function of time, calculated by Thompson et al. [IB ]. The 
progenitor mass is 15M Q . 

neutrinos and antineutrinos. 

The neutrino luminosity curve is quite characteristic 
among the different flavors. In particular, there is a very 
sharp peak of v e called the neutronization burst, whose 
duration is typically ~ 10 ms and peak luminosity is 
~ 10 53 erg s _1 . This strong peak is illustrated as follows. 
As a supernova shock moves outward, it dissociates nu- 
clei into free nucleons, which triggers the deleptonization 
process e~p — > v e n; these v e build up a sea because they 
are trapped and advected with matter. When the shock 
crosses the v e neutrinosphere, within which the created 
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FIG. 2: Original neutrino spectrum integrated to 0.25 s after 
core bounce, calculated by Thompson et al. |T3| . The progen- 
itor mass is 15M©. 



v e are trapped, they are abruptly emitted. For the other 
flavors v e and v x , there is no such sudden burst; both 
luminosities glow rather gradually and they are similar 
to each other. 

The other characteristic that provides information on 
the flavor conversion mechanism as well as neutrino decay 
is the hierarchy of the average energy (E Ue ) < (-Ep c ) < 
(_E„ x ) as clearly seen from Fig. [21 neutrino flavor conver- 
sion and decay also change the spectral shape. This en- 
ergy hierarchy is explained as follows. Since v x interacts 
with matter only through neutral-current interactions in 
supernovae, they are weakly coupled with matter com- 
pared to v e and v e . Thus the neutrinosphere of v x is 
located deeper in the core than that of v e and v e , which 
leads to higher temperatures for v x . The difference be- 
tween v e and v e comes from the fact that the core is neu- 
tron rich and the v e couples with matter more strongly, 
through the v e n — > e~p reaction. 

Although we only use one specific model the sig- 
nature of the neutronization burst appears in all reliable 
numerical simulations. It is quite natural that the height 
and width of such a peak are dependent on the super- 
nova parameters as well as numerical approaches. How- 
ever, we are not interested in such a slight difference but 
focus only on the rough signature — i.e., the appearance 
or absence of a neutronization peak in the v e channel. 
We cannot discuss in any details without a concrete and 
reliable original model, but even such rough information, 
if ever detected, would bring a very fruitful and novel 
perspective to particle physics. As for the neutrino spec- 
trum, there is a general tendency of the hierarchy of the 
average energies — i.e., (E Ve ) < (Ep e ) < {E Vx ) as already 



noted and this tendency seems to be common to almost 
all numerical calculations. However, the spectral shape 
and especially the average energy ratio between v e and 
v x are still matter of controversy, and we cannot conclude 
that the spectral information at the detectors would be 
very useful at present (see , for a comparison among sev- 
eral calculations, Ref. [13] ). Contrary to the luminosity 
curve, the neutrino spectrum does not indicate such a 
prominent signature as a neutronization burst, and there- 
fore, the obtained spectrum would be of secondary im- 
portance for the purpose of this study. In the near future, 
however, the situation may become significantly better, 
especially if some process relevant for successful super- 
nova explosions were discovered in computers; in that 
case, the spectrum as well as luminosity curve would be 
very useful to precisely obtain the unknown property of 
neutrinos. 

Before moving on to a discussion including neutrino 
decay, here we shortly describe the flavor conversions in- 
side the supernova envelope without any other processes 
such as decay. Neutrinos of different mass eigenstates 
are expected to be incoherent with each other when they 
reach the detector, and then the number intensity (i.e., 
number per area, time, energy, and solid angle) of 9 e can 
be simply represented by 

I» e (L) = |t/ e i| 2 / Pl (i?SN) + |f/e 2 | 2 ^ 2 (i?SN) 
+ \U e3 \ 2 h 3 (R SN ) 
~ cos 2 6>i 2 / Pl (i? SN ) +sm 2 e»i 2 / P2 (i? SN ), (6) 

where L and i?sN represent the distance to and radius of 
the supernova, respectively. The second equality comes 
from the fact that the value of #13 is constrained to be 
negligibly small from reactor neutrino experiments [29^ 
and #23 — tt/4; from the solar and reactor neutrino obser- 
vation the obtained value for #12 is cos 2 612 ~ 0.7 (LMA 
solution) 0, • The intensity at the supernova surface 
Iv{Rsn) reflects the flavor conversions during propaga- 
tion inside the supernova envelope. Flavor conversions 
during neutrino propagation have been extensively stud- 
ied by many researchers (see, e.g., Ref. |3(j), but we 
briefly summarize them here. In the case of a normal 
mass hierarchy {mi <^ 7723) , on which we focus in this 
study, the intensity of each mass eigenstate at the sur- 
face is fairly well known to be 

IpARsn) = -Mo), 

h 2 (R S N) = 7 C3 (i? SN ) = ^(0), (7) 

where represents the neutrino intensity of each fla- 

vor eigenstate at the neutrinosphere, for which we use 
the results of numerical simulation shown in Figs, ^and 
El By a combination of Eqs. JHl and (0), we can see that 
owing to flavor conversion inside the supernova, about 
30% of the detected v e would originate from v x at pro- 
duction; this would harden the obtained spectrum at the 
detectors. On the other hand, for the neutrino sector, 
the final expression for the intensity of v e is given by the 
same expression as Eq. (JHJ but with replacing i> by v. 



5 



Flavor conversions inside the supernova envelope are, this 
time, a little bit complicated; the unknown parameter 
6*i3 strongly affects the results. Instead of #13, we rather 
use so-called "flip probability" at the higher resonance 
point P H [13, which equals 0(f) when sin 2 26>i 3 > fCT 3 
(sin 2 2#i3 < 10 -5 ). The expressions corresponding to Eq. 
(|7|) for the neutrino sector are then given by 

M^sn) = /„.(0), 

UARsn) = Pjr^.(0) + (l-P ff )J I/ .(0), 
L 3 {Rsn) = {1-P h )UM+PhLA0). (8) 

These expressions are necessary for estimating the neu- 
trino flux in the case of possible decay, since 14 — > 9j 
might occur. 

IV. FORMULATION 

In this section, we derive new formulation for the de- 
tected v e flux, which includes both flavor conversion and 
decay during propagation. The v e intensity at the detec- 
tors is represented by 

I Pe (L,E) = \U el \ 2 I Pl (L,E) + \U e2 \ 2 h 2 (L,E) 



+ \U e3 \ 2 I* 3 (L,E), (9) 

which is similar to Eq. ©, but the intensity of the spe- 
cific mass eigenstate is no longer conserved during its 
propagation owing to decay, Ip t (L) ^ I„ i (-Rsn)- Here and 
from this point on, we explicitly show the neutrino energy 
E. The effect of neutrino decay on the intensity of each 
mass eigenstate T>i is included by adding the appearance 
and disappearance terms to the transfer equation — i.e., 

= (r,2J) 

j<i 

+ V / dE' [ME', E)r, ] ^(E')I P] (r, E') 

+iP f (E',E)r Uj ^ i (E')I Uj (r,E')] , (10) 

where we define Jjj = r c , i ^ tUj + etc. A similar 

formulation holds for the neutrino sector, although we 
do not show it explicitly. The first and second sums of 
Eq. I|10[) reflect the disappearance and appearance of 
&i, respectively. Fortunately, this set of formulations can 
be analytically integrated from i?gN to L. In the case 
of three-flavor context and normal mass hierarchy, the 
solution to Eq. (fTUfl is given by 



hAL,E) = 



I Pl (R SN ,E)+ / dE' 
Je 



_ e -(r 13 (E')+r 23 (E'))L 

r 13 (E>) + r 23 (E>) 



-J 3 ^(E',E) 



1 _ e -r 12 (E')L 

r 12 {E>) 



J 2 ^i(E',E) 



dE" 



-r 12 (E')L 



h 2 (L.E) 
h 3 {L.E) 



r 12 {E' 
e- r ^ L h 2 {R S n,E) + 



dE' 



_ e -(r 13 (E")+r 23 (E"))L 

r 13 (E") + r 23 (E") 



e -r 12 (E)L _ e -(r 13 (E')+r 23 (E'))L 



r 13 {E') + r 23 (E') ~ r 12 {E) 



Jz^2^i (E , E' , E) 
J 3 ^ 2 (E',E), 



e -(r 13 (E) + r 23{ E))L Ip3(Rs ^ E)i 



where 



J3^2^i{E" , E' , E) — 



1 



r l3 {E") + r 23 {E") - r 12 (E') 



x ME', E)r 92 ^ Pl (E')J 3 ^ 2 (E", E') + iP f (E> \E)r V3 ^ Pl (E')J 3 ^ 2 (E", E') 
J^E'.E) = ME',E)r^ Pj {E')h z {R SN ,E')+^ f (E' 1 E)r^ j (E')I^(R SN ,E'), 
J^(E', E) = ME', E)r v% ^ V] {E')I Vi (Rsn, E') + i> f (E', E)E- V ^ V] (E')h z (R SN , E'). 

I 



(11) 

(12) 
(13) 



(14) 
(15) 
(16) 



In these expressions, we have used the assumption that 
EijRsN •C 1; i.e., neutrinos never decay during their 
propagation inside the supernova envelope. With this as- 
sumption, the intensity at the stellar surface I„(Rsn, E) 
is, also in this case, represented by Eqs. Q and JSJ. 
Thus we obtain the intensity of v e at the detector using 
Eqs. ©, dHJ)— <dSJ) , 0, and ©. 



Although we have given a quite general expression, we 
are rather interested in the more specific case, in which 
the neutronization peak appears in the v e channel; this 
is realized when the neutrino masses are strongly hier- 
archical as already discussed in Sees. Inland IIIII In this 
case, we obtain r Pi ^ Uj — r^—^. = -Tji/2, etc., and the 
energy distribution factions are given by Eq. Q. From 
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this point on, we use Tw/m, Tis/m, and Ti^/m as free pa- 
rameters, which are related to JJy by Eij(E) = m/fajE). 



RESULTS 



The obtained number flux and fluence (time- integrated 
flux) of D e , using Eqs. ©, dJ-lJISl, 0, and © are 
shown in Figs. I3HE1 In the inset of each figure, we show 
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FIG. 3: Energy-integrated flux at the Earth in the case of 
adiabatic conversion (Ph = 0), for various values of nz/m. 
The other parameters (riz/m and T2i/m) are set to infinity 
The expected number of events at the detector of 640 kton is 
shown in the inset. 



the expected number of detection at the water Cerenkov 
detectors with a fiducial volume of 640 kton [20 times 
larger than that of the Super-Kamiokande (SK) detec- 
tor], by assuming that the supernova occurred at 10 kpc. 
The values labeled in the vertical axis could be rescaled 
by using a factor of (10 kpc/Z?) 2 (Ffid/640 kton) in other 
cases, where D represents the distance and Vad the fidu- 
cial volume. The cross section for the dominant catching 
process D e p — > e + n is fairly well understood and we used 
that given in Ref . jU . In addition we have used the trig- 
ger threshold expected to be installed to SK-III (after full 
repair) ; with this threshold the electrons and positrons of 
the energy of 3 MeV can be detected at 100% efficiency. 
We neglect all other processes such as electron scattering 
because of their subdominance. 

Figures |3| and 0] indicate the energy- and time- 
integrated flux, respectively, in the case of the adiabatic 
conversion, i.e., Ph = 0. The fluxes are evaluated for 
various values of Tiz/m with t^/to and Ti^fm fixed to 
infinity. The shape of the luminosity curve is found to 
strongly depend on the lifetime of the 1^3 (D3) — > T>\ mode. 
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FIG. 4: Fluence (time-integrated flux) at the Earth in the 
case of adiabatic conversion (Ph = 0), for various values of 
m/m. The other parameters ("m/m and T23/7T1) are set to 
infinity. The expected number of events at the detector of 
640 kton is shown in the inset. 



This is because in the case of adiabatic conversion, the 
original v e becomes 1/3 at the stellar surface [Eq. (JSJ] and 
they decay into D\, which dominantly couple to v e . Thus, 
the peak of the neutronization burst clearly appears at 
the detectors for ri 3 /m < 10 5 s/eV. As for the spec- 
trum, the energies of daughter neutrinos are significantly 
degraded as shown in Fig. 0] and give a very characteris- 
tic signature. However, since the cross section is roughly 
proportional to E 2 and highly insensitive to low-energy 
neutrinos, the expected event number is almost the same 
at such a low-energy region as shown in the inset of Fig. 
01 This degradation of the neutrino energy due to its de- 
cay also causes the actually detected neutronization peak 
to be less prominent compared to that seen in the flux. 

The case of nonadiabatic conversion (Ph = 1) is shown 
correspondingly in Figs. and This time, the rele- 
vant parameter is changed to Tn/m because v e created 
by the neutronization is converted into at the stellar 
surface. The characteristics appearing in both the lu- 
minosity curve and spectrum are similar to those in the 
case of adiabatic conversion, but the total event number 
is slightly smaller. This difference comes from the fact 
that the detected v e is coupled to by ~ 30% [Eq. 
and the disappears owing to its decay. On the other 
hand, the disappearance of D3 does not directly affect the 
expected event number since they hardly couple to v e . 

In order to discuss the parameter dependence of this 
mechanism, we simply define the following quantity, 
which represents the peak sharpness of the time profile. 
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FIG. 5: The same as Fig. EJbut for the case of nonadiabatic 
conversion (Ph = 1), for various values of T12/TO with fixed 
T13/771 and T23/m to infinity. 



the neutronization burst is more prominent. We plot the 
contour map of Rt against the values of r 12 /m and r 13 /m 
assuming several values for T23/m; the result in the case 
of adiabatic (nonadiabatic) conversion is shown in Fig. [7| 
(Fig. |8J. In both figures, the solid, dashed, and dotted 
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FIG. 6: The same as Fig. [I] but for the case of nonadiabatic 
conversion (Ph = 1), for various values of i~i2/m with fixed 
T13/771 and T23/m to infinity. 



Namely, it is defined as 

event number in region (I) 



Rt 



event number in region (II) ' 



(17) 



where regions (I) and (II) are defined in the insets of Figs. 
|3and[SJ A larger value for Rt means that the peak of 



FIG. 7: Contour plot of Rt against the (riz/m, Tis/m) plane, 
in the case of the adiabatic conversion (Ph = 0). The val- 
ues of T23/TO are taken to be 10 7 (solid curves), 10 3 (dashed 
curves), and 1 (dotted curves) s/eV. Each curve of the same 
type is equally spaced by the value of Rt with the indicated 
largest and smallest values. 



curves represent the case that the values of T23/m are 
10 7 (essentially no decay due to the 3^2 mode), 10 3 , 
and 1 s/eV, respectively. Each curve of the same type 
is plotted at an equally spaced level, while the largest 
and smallest values are indicated. It can be clearly seen 
that in the case of Pjj — 0, the Rt is strongly dependent 
on the parameter T13/771 but highly insensitive to Txz/m. 
This point also holds for the nonadiabatic case with the 
corresponding replacement, T13 <-» T12. 



VI. DISCUSSION 

In order to discuss how to discriminate one from the 
other decaying models, in addition to Rt, we use another 
quantity that represents the spectral hardness: 



Re 



event number for E e > 25 MeV 
event number for E e < 15 MeV 



(18) 



The values for Re are obtained from the detected 
spectrum — i.e., from the insets of Figs. 0] and El As 
already discussed in Sec. IIIII because the average energy 
difference between each flavor neutrino is still a matter 
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FIG. 8: The same as Fig. 0but for the case of the nonadia- 
batic conversion (Pr = 1). 



of controversy, the obtained spectrum would be also af- 
fected by such uncertainties. We believe, however, that 
our treatment is quite reasonable, since using the simple 
quantity Re would make the discussion rather insensitive 
to such unsettled details. 

Figure |5] shows the location of each model on the 
(Re,Rt) plane for the adiabatic case. We also show 
the la statistical error bars of Re and Rt', the size of 
these errors changes as {D/10 kpc)(Vfid/640 kton) -1 / 2 
in other cases. Labels such as "2 — > 1" represent the 
relevant decaying mode, while the other modes are as- 
sumed to be stable; the label "All" represents the case 
that T12 /m = Ti3 /m = T23 / m. Points of the same symbol 
show how their location changes with lifetime; each point 
of one symbol represents a model with a specific value of 
lifetime-to-mass ratio r/m, which is equally spaced (log- 
arithmically) from 10 4 to 10 7 s/eV. The dotted lines con- 
nect these points just to guide our eyes. All the modes 
are degenerate when r/m — 10 7 s/eV at Re = 0.23 and 
Rt = 0.27, which means that there occurs essentially no 
decay. From this figure, models with extreme parameter 
values can be distinguished from one another, if the cur- 
rently planned megaton-class detectors, such as Hyper- 
Kamiokande and UNO, detected the galactic supernova 
neutrino burst. On the other hand, in the case of cur- 
rently working detectors such as SK, the errors become 
very large by a factor of > v20, and even using these very 
simple quantities Re and Rt, it would be quite difficult 
to derive some information. 

We here briefly illustrate the behavior of each model, 
shown in Fig. [5] As already discussed in the previous 
section, the decaying mode from ^3(1^3) to Di makes the 
value of Rt larger owing to the appearance of a neu- 



FIG. 9: Neutrino decay model plotted on the (Re,Rt) 
plane in the case of adiabatic conversion. The error 
bars include only statistical errors, and are at the la 
level, but their size should be accordingly rescaled as 
(D/10 kpc)(Vfi d /640 kton)" 1/2 . The labels represent the rel- 
evant mode (see text) and points of the same symbol indicate 
the model with a specific value of r/m, which is equally spaced 
from 10 4 to 10 7 s/eV, logarithmically. 



tronization peak. A similar explanation applies to the 
3 — > 2 mode but its prominence is reduced because the D2 
state included in the u e is smaller compared to v\ state. 
The 2 — ► 1 mode does not change the peak sharpness 
Rt, since in the case of adiabatic conversion, the 1*2(^2) 
at the stellar surface does not contain any component 
from the original v e . Instead of an almost constant Rt, 
the spectral hardness Re significantly changes with the 
value of T\ijm. This is also easily illustrated as follows. 
At the supernova surface, the vi and Vi both originate 
from v x , which shows the hardest spectral shape. With- 
out any decay, the obtained v e signal then contains an 
~ 30% amount of the original v x . On the contrary, if 
the decaying mode 2 — > 1 were relevant, the original v x 
component, which otherwise should contribute to the v e 
spectrum, would disappear and instead it would appear 
as v\ but with a significantly reduced energy; this makes 
the spectral hardness Re considerably small. The "All" 
model, in which we assumed 712/m — Tia/m = T23/m, 
includes both effects given above; i.e., the peak sharpness 
Rt increases owing to the 3 — ► 1 and 3^2 modes, while 
the spectral hardness Re decreases owing to the 2 — > 1 
mode. 

Figure ^3 is the same as Fig. but for the case of 
nonadiabatic conversion Ph = 1. In this case, the de- 
caying mode from ^3 and D3 does not essentially change 
the expected signal, because they are not coincident with 
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decaying models for the adiabatic conversion with rjvn — 10 4 
s/eV are also plotted for comparison. 



the original v e at the stellar surface as well as having es- 
sentially no coupling to v e . The original v e , instead, in 
this case, appears as vi\ thus the decaying mode 2 — > 1 
considerably changes the detected signals. 

There also exists another mechanism that possibly 
changes the original v e into a detected D e , resulting in 
the appearance of a sharp peak due to a neutronization 
burst at the detectors — i.e., resonant spin- flavor (RSF) 
conversion [TTL fl^ . This mechanism is induced by the 
interaction between a supernova magnetic field and the 
Majorana magnetic moment of neutrinos. According to 
Ref. [Tl|, the very sharp peak of a neutronization burst 
could appear owing to the combination effect of the RSF 
and ordinary matter-induced conversion, if the following 
conditions are all satisfied: (i) the mass hierarchy is in- 
verted, (ii) the value of #13 is sufficiently large, and (hi) 
the neutrino magnetic moment as well as supernova mag- 
netic field is large enough to induce the adiabatic RSF 
conversion (but see also Ref. H2)- In order to compare 
the RSF mechanism with the decaying models, we plot 
the model groups given in Ref. ^lj m Fig- E] F° r com- 
parison we also plot the decaying models with r jra = 10 4 
s/eV and Pjj = 0. The models with r/m = 10 7 s/eV (no- 
decay model) are degenerate with group A due to the 
RSF conversion. This figure clearly indicates that the 
RSF mechanism potentially gives a far more character- 
istic T> e signal at the detectors; both the peak sharpness 
of time profile Rt (group C) and spectral hardness Re 
(group B) are very prominent, and even the SK detector 
enables us to discriminate these model groups. 



VII. CONCLUSIONS 

Nonradiative decay of neutrinos is not constrained suf- 
ficiently; the most stringent lower limit to the lifetime t\i 
is obtained from the solar neutrino observation but it is 
very weak 0. Thus, neutrino decay possibly affects the 
neutrino signal from other astrophysical objects such as 
supernovae. 

Using newly derived formulation, which includes both 
flavor conversions inside the supernova envelope and neu- 
trino decays during propagation in vacuum, we calcu- 
lated the expected neutrino luminosity curve as well as 
the spectrum at future large volume water Cerenkov de- 
tectors. In these calculations, we particularly focused on 
the decaying model such that the original v e appears in 
the v e signal as the result of flavor conversion and decay. 
This is because such a situation can give the appearance 
of a sharp peak in the time profile due to the neutroniza- 
tion burst, and it could be easily recognized. We discuss 
that this actually may be realized if the neutrino masses 
are strongly hierarchical, and we have assumed it in ac- 
tual calculations. The lifetimes of three decaying modes, 
Ti2) r i3> an d T23, are taken to be free parameters, and 
the cases of adiabatic and nonadiabatic conversions are 
treated independently. The results of the calculations 
are shown in Figs. I^HHland the neutronization peak can 
be significantly prominent in future megaton-class water 
Cerenkov detectors. 

In order to discuss the parameter dependence of the 
neutrino signal, we introduce the rather simple quanti- 
ties Rt and Re , which represent the peak sharpness and 
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spectral hardness, respectively. As shown in Figs. [7] and 
|SJ the value of R T is strongly dependent on the value of 
P/f as well as the relevant lifetime. From Figs. l^land lTUl 
we see that the behaviors of each model on the (Re, Rt) 
plane are considerably different from one another. But 
we also show that the location of these decay models 
on this plane does not give prominent properties enough 
for us to distinguish using current detectors such as SK, 
on the contrary to a very significant dispersion due to 
the RSF conversion (Fig. Ill|) . Finally, we again stress 



that the appearance of the neutronization peak clearly 
indicates nonstandard properties of neutrinos; neutrino 
decay would then be one of the possible candidates. 
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